Kinetic equation approach to diffusive superconducting hybrid devices 
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We present calculations of the temperature-dependent electrostatic and chemical potential dis- 
tributions in disordered normal metal-superconductor structures. We show that they differ appre- 
ciably in the presence of a superconducting terminal and propose an experiment to measure these 
two different potential distributions. We also compute the resistance change in these structures due 
to a recently proposed mechanism which causes a finite effect at zero temperature. The relative 
resistance change due to this effect is of the order of the interaction parameter in the normal metal. 
Finally a detailed calculation of the resistance change due to the temperature dependence of An- 
dreev reflection in diffusive systems is presented. We find that the maximal magnitude due to this 
thermal effect is in general much larger than the magnitude of the novel effect. 

PACS numbers: 74.50. +r;74.80.Fp 



I. INTRODUCTION 

Mesoscopic structures in which normal metal wires or 
semiconductors are attached to superconductors have re- 
ceived a fair amount of attention in the past few years. 
In particular devices known as Andreev interferometers, 
in which two superconducting terminals with different 
macroscopic phases are present, have been in the focus 
of interest. The conductance of these structures, in which 
electrons and holes undergo multiple Andreev reflection 
depends on the phase difference of the connected super- 
conductors, hence the name Andreev interfero*netry. 

Since the prediction of Andreev reflection,!!! the the- 
ory of charge transport through N-S junctions has been 
well established. However, the practical implications 
of this phenomenon for the sophisticated nano-structures 
that can nowadays be realized are not always clear. The 
reason for this is the coherent nature of multiple Andreev 
reflection which determines the physical behavior of these 
devices. These new technological developments resulted 
in the current revival of the topic in mesoscopic physics. 

In the last few years a large number of Andreev-. ip=. 
terferometers have he <p, studied both theoreticallyETtEl 
and exp eriment allyO~liia Particularly the experiment of 
Rcf. |l4] motivated the research presented here. In this 
experiment, the resistance of a cross shaped diffusive nor- 
mal metal was measured. The two branches of the cross 
perpendicular to the current path were in contact with 
a large superconducting loop. The phase difference be- 
tween the superconducting end points of the loop could 
be controlled by a small current through the loop or, al- 
ternatively, by applying a magnetic field. The resistance 
of the structure oscillated non-harmonically as a function 
of the phase difference by about 10% of the normal state 
resistance. These results were unexpected because in 
the conventional theory of the proximity effect, in which 
electron-electron interactions in the normal metal region 
are disregarded, the zero-voltage, zero-temperature re- 
sistance of a diffusive metal is predicted to be phase- 



independent .Il§li3 Furthermore, the large amplitude and 
the observed 2tt periodicity ruled out the possibility of a 
weak localization effect, since the latter was predicted to 
shaw-a phase-dependence with a fundamental period of 
7r.EZrEa To this day, resistance oscillations with 7r period- 
icity remain unobserved. 

Recently we proposed a new mechanism which pro- 
vides a phase-dependent_resistivity in a diffusive conduc- 
tor at zero temperature. lj This scheme takes into account 
the fact that the electron-electron interaction induces a 
weak pair potential in the normal metal. As a result, 
Andreev reflection occurs in the entire structure, rather 
than only at the N-S interfaces. This results in a phase- 
dependent resistance change which is proportional to the 
interaction parameter A in the normal metal and can be 
of either sign, depending on the sign of A. Although the 
experiment of Ref. ^ could be explained in terms of the 
proximity effect theory and the results were shown to 
be caused by the finite .temperature at which the exper- 
iments were pcrformedjlj it would be challenging to ob- 
serve the resistance oscillations predicted in Ref. [H| This 
would also be of practical interest since it would provide 
the means to directly measure the interaction parame- 
ter in the normal metal. However, besides the fact that 
electron-hole coherence influences the resistance, it also 
manifests itself in a non-trivial distribution of the elec- 
trostatic and chemical potentials in the structure as we 
will show below. 

The remainder of this paper is organized as follows. In 
section || we briefly discuss the influence of phase coher- 
ence on transport properties and potential distributions 
in small diffusive structures. Section III contains the the- 
oretical foundation of our calculations. We first review 
the relevant techniques of the Keldysh formalism for dif- 
fusive superconductors and then derive the equations for 
the Green functions and distribution functions which de- 
termine the electric transport properties of the system. 
The next three sections are devoted to several applica- 
tions of the theory. We first calculate the temperature- 
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dependent electrostatic and chemical potential distribu- 
tions in a simple one dimensional structure in section IV . 
A second application is presented in section [v|, where we 
calculate the resistance change at zero temperature due 
to the induced pair potential in the normal metal region 
for two experimentally relevant geometries. The third 



and last application is discussed in section VI. There we 
focus our attention on the influence of a finite tempera- 
ture on the resistance in these structures. Some of the 
results in sections M and VI were published in a prelimi- 
nary form in Rcf. [TBI However, here we additionally give 
a detailed description of the performed calculations. We 
summarize our conclusions in section VII. 



a) 
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II. COHERENCE EFFECTS IN ULTRASMALL 
DISORDERED STRUCTURES 

Owing to the advance in nanofabrication techniques 
over the past years, the fabrication of hybrid metallic su- 
perconducting structures with a characteristic size of a 
few microns or less has nowadays become possible. If 
these small structures are at a sufficiently low tempera- 
ture, the quasiparticles in the metal can no longer pen- 
etrate into the superconductor due to the large super- 
conducting gap. As a result, the lowest order process 
that determines the resistance of the system is Andreev 
reflection, in which an electron is reflected as a hole or, 
alternatively, in which an electron pair enters the super- 
conductor. This reflection causes electrons and holes in 
the diffusive metal to be p hase coherent over distances 
of the order of £ = \JT>/T ^> I, where T> = ^v-pl is the 
diffusion constant and I is the elastic mean free path. 

This phase coherence between electrons and holes dras- 
tically alters the physics of transport through such sys- 
tems. The most striking feature is that the electrostatic 
potential and nonequilibrium chemical potential are no 
longer distributed linearly through the sample. The non- 
linearity of the electrostatic potential implies a nonuni- 
form resistivity distribution and consequently a non-local 
resistance of the structure. This nonlocality is a fun- 
damental feature of the coherent nature of Andreev re- 
flection. Moreover, at finite temperatures the transport 
properties of the system cease to be distributed uniformly 
over all energies. Hence, a calculation of these quantities, 
to which the main part of this article will be devoted, 
must first consider them at each energy individually and 
then integrate over all energies. 

Another manifestation of the phase coherence in the 
normal metal is the difference in the distribution func- 
tions of the electrostatic and chemical potentials. In a 
normal system both would be equal, but this changes 
when one of the leads is brought into the superconduct- 
ing state. Whereas the former is simply determined by 
the distribution of charge in the system, the latter can 
only be defined for small perturbations from equilibrium, 
i.e. when the quasiparticle energies are much smaller 



b) 

FIG. 1. Schematic setup of an experiment to measure the 
electrostatic and chemical potential: (a) Top view, (b) Cross 
section along the dashed line in a). For details see text. 



than the superconducting gap. To show how this defini- 
tion comes about we consider the normal current through 
a disordered N-I-S junction. Zhou, Spivak, and Zyuzin 
showed that for small quasiparticle energies e this current 
can be written in the following wayO 

j(e,x N ) = t{f T (e,x s ) - f T (e,x N )} F(e, x N ,x$), (1) 

where t is the transparency of the tunnel barrier, xn and 
Xg denote the normal-metal side and superconducting 
side of the barrier, F is some function of e, in and xg, 
and /t is the nonequilibrium distribution function. This 
equation shows that at low temperatures the nonequilib- 
rium chemical potential can be associated with the distri- 
bution function and consequently that it is a measurable 
quantity. As will be discussed in section [V, the electro- 
static potential decreases faster than linear in the vicin- 
ity of the superconductor due to the decreased density 
of states near the N-S interface. In contrast, the chem- 
ical potential changes only a little in the presence of a 
superconducting terminal and consequently the ratio of 
the electrostatic and chemical potential vanishes near the 
superconductor. 

A possible experimental setup to measure the differ- 
ence between the electrostatic and chemical potential is 
drawn in Fig. |l|. The figure shows a diffusive wire con- 
necting two reservoirs (not shown) in the presence of two 
different probes. The voltage probe measures the chemi- 
cal and the capacitance probe the electrostatic potential. 
As shown in Fig. |l](b), the voltage probe consists simply 
of a metallic lead separated from the wire by a thin oxide 
layer. The latter is indicated by the dark shaded region. 
The capacitance probe is slightly more complicated. In 
principle one could use a single metallic gate separated 
by a thick oxide layer to reduce tunneling from lead to 
wire. However, because such a gate would not only cou- 
ple capacatively to the charge in the wire directly un- 
derneath the tip but also strongly to the surroundings, 
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we propose a slightly different method: A small metallic 
island is deposited on the edge of the wire as indicated 
in the figure. This island is weakly coupled to two extra 
leads through which a current can flow. In this case the 
electrostatic potential capacatively induces a charge on 
the island, thus very sensitively changing the measured 
Coulomb threshold. Both probes can be calibrated since 
they should measure the same potentials if both reser- 
voirs are in the normal state. If one of the reservoirs 
is brought into the superconducting state, both poten- 
tials should change. If one would attach different probes 
along the wire it would be possible to measure the poten- 
tial landscapes in the wire. In section IV we will return 



to this subject in a more quantitative fashion, but we will 
now first discuss the necessary theory. 



III. METHOD 

A. Kinetic equations for the distribution functions 
in the dirty limit 

To describe the system we use the nonequilihrium 
Green function method first derived by KeldyshEJ and 
later further developed for superconductivity by Larkin 
and Ovchinnikov.EJ Although this framework is rather 
formal it has one big advantage over e.g. the scattering 
approach. As was mentioned above the resistance of the 
structure is nonlocal. However, using the Keldysh tech- 
nique it is possible to express scattering processes in the 
structure as well as other relevant physical quantities, 
like the resistance, in terms of local Green functions in 
coincident points. This property of the formalism simpli- 
fies the calculations considerably. To establish notation 
and to remind the reader of the basic theory we briefly 
review the essential ingredients of the Keldysh formalism 
and the quasiclassical approximation for diffusive super- 
conductors. For more extensive reviews we refer to e.g. 
Rcfs. Eland |23. 



In this formalism the Green function is given by the 
(4 x 4) matrix 



G = 



G A G K 
G R 



(2) 



where G A ,G R , and G K are the advanced, retarded, and 
Keldysh Green function, which are (2 x 2) matrices in 
Nambu space given bycJ 



G A (l,f) 



G A (1,1') F A (1,1') 
F A t(i,i') -G A (1',1) 



and analogous equations for G R and G K . Throughout 
this paper the symbol "check" will be used to denote 
(4 x 4) matrices and "hat" for (2 x 2) matrices. The 
normal and anomalous Green function are given by: 



G R (1, 1') = -%e{h - fiO#(l),^t(i')] + ), (4a) 

G A (1, 1') = iO(t v - fi)<[V(l), V f (!')]+>, (4b) 

G K (l,lO = -i(^(l),V t (lO]) I (4c) 

F R (1, 1') = -i6{h - ty)([^(l)^(l')} + ), (4d) 

F A (l,l / )=ifl(ti'-*i)<[V(l),V'(l / )]+> ) (4e) 
F K (l ! lO=i([^(l),^(l')]), (4f) 

where ip(l) — i/j(ti, ri) is the electron field operator. 

We proceed by introducing the center-of-mass and rel- 
ative coordinates r = |(ri +rj/) and r' = ri — ry and by 
Fourier transforming the Green function with respect to 
the relative coordinate: G(r, p) = J dr' exp(— ipr') G(r+ 
|r', r — |r'). We apply the quasiclassical approximation, 
which is based on the fact that the Fermi energy in the 
system is much larger than all other energy scales. This 
means that all relevant physical quantities vary spatially 
on a length scale that is much larger than the Fermi wave- 
length. In this case it is useful to introduce the so-called 
quasiclassical Green function g which is integrated over 
£ — — 



ft: 



g(r,p,ti,t v ) = - / tZ£ p G(r,p, 



(5) 



Here p in the left hand side denotes the fact that the mo- 
mentum dependence of the quasiclassical Green function 
is restricted to dependence on the direction of p only. In 
this approximation the magnitude of the momentum is 
fixed at |p| — pp. This quasiclasissical Green function 
satisfies the normalization condition 



dti»ff(*i,ti»)fl(*i",ti') = -*!')■ 



(6) 



In the case of a superconductor with short elastic mean 
free path, i.e. in the diffusive regime, it is feasible to 
expand thp-,Green function to first order in spherical 
harmonics :Ea 



g = g s + pg P , pgp^g s 



(7) 



where the functions g s and g p no longer depend on the 
direction of p. Using the normalization condition (^|) we 
find an expression for g p , which is then substituted back 
into Eq. (Q). The thus obtained Green function is then 
averaged over all angles of p. In the stationary case, the 
Green function depends on the time difference r = t\ —ty 
only. Performing the Fourier transform with respect to 
this time difference, we obtain the equation of motion for 
the Green function g £ c?rg s (r) exp(?£r) (we drop the 
' n)E3 



(3) index "s" from now on)c 



■ 2?V (ge^ge) + i[H, 9e] + g e ] = 0, 

I, 



9l 



(8a) 
(8b) 



where T> = ^vfI is the diffusivity, / is the elastic mean 
free path, H = eipl + ea z — A and A and a z are given 
by: 
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A 
A 



A 

-A* 



where A is the pair potential in the metal and 





( a z 









*. = l 


v o 


Oz J 







(9) 



(10) 



In Eq. ( [8a| ) elastic non spin-flip impurity scattering has 
been taken into account in the Born approximation, caus- 
ing the presence of-,the elastic mean free path / in the 
diffusion constants Hence, the self energy matrix £ in 
Eq. ([Sal), which has the same structure as the Green 
function (^), takes into account processes like spin- flip 
scattering and the inelastic scattering of electrons with 
phonons and (magnetic) impurities. 

The general expression for the electrostatic potential 
ip in Eq. JSaL), which follows from electroneutrality in the 
metal, isEjxJ 



tp(x) 



1 



de Tr gf(x) 



(11) 



Throughout this paper the electrostatic potential is as- 
sumed to be time- independent. 

The advanced and retarded Green functions determine 
the dispersion of the quasiparticles. However, to solve 
a transport problem we need to know how the energy 
spectrum is filled by extra quasiparticles when the system 
is driven out of equilibrium. This is determined by the 
Keldysh component g K of g, which can be expressed in 
the advaneed-and retarded ones using two distribution 
functions t2rE3 



/ = /l1 + Jw*. 



(12a) 
(12b) 



In a spatially slowly varying electro-magnetic field, the 
equations for the two distribution functions are (dropping 
collision integrals because they account for inelastic scat- 
tering processes, and time derivatives because we seek to 
find stationary solutions only) 



(13a) 



(13b) 



PVTr {V/ L (l - gfg^)} + VW (f Tje ) = 
PVTr {V/ T (l - a z gfa z g^)} + VVf h j e + 
2 i / T Tr{a R +3 £ A ) A] 

where j e = Tra z {gfdgf - g £ dg^}. 

To close the set of equations we finally need an equa- 
tion for the pair potential in the normal metal region. 
This expression for A can be derived from the selfconsis- 
tency relation: 

Mr) = ^Jds {g?(v,p,r)} od , 

a R -.9e A } .d> ( 14 ) 



Ti 



de tanh 



where A = gN(0) is the interaction parameter, g, times 
the density of states at the Fermi level, N(Q). The sub- 
script o.d. denotes the off-diagonal part. 

This concludes the derivation of the distribution func- 
tions. We now have a closed system of equations that 
in principle must be solved sclfconsistently. In the next 
section we will discuss specific circumstances that allow 
for a simplification of the equations, enabling us to solve 
them perturbatively. 



B. Approximations 

In this section we discuss the assumptions and approx- 
imations that enable us to simplify the theory. We sub- 
sequently derive the final set of equations that we will 
use. 

We start by noting that Eq. J8a| ) for the Green func- 
tion still contains the self energy matrix S. However, be- 
cause we are only interested in the case where the phase 
breaking length is much larger than the system size, it is 
reasonable to disregard inelastic scattering processes and 
hence, we neglect X from now on. Eq. ( |8a| ) then reduces 

tcM ' 



VV (g £ Vg £ )-i[H,g £ }=0. 



(15) 



We now have equations for the diagonal components g^ 
and g^ of the Green function (Q). 

We parametrize the advanced Green function in the 
following way: 



cos f 



ie 1 ^ sin( 
— cos# 



(16) 



thus ensuring that g £ = 1. In general cj> and 9 are com- 
plex and depend on energy and position. In a structure 
with two superconducting terminals, g^ will depend on 
the phase difference <\>\ — <f>% between the two supercon- 
ductors. However, if only one superconducting reservoir 
is present, the resistance of the structure will not depend 
on the absolute phase and we can put tf> = 0. 

In the case of sufficiently small quasiparticle, thermal 
and Thouless energies; e,kBT,T>/L 2 -C As, where As 
is the energy gap in the superconductor, the advanced 
Green function can be written in the following way:Ell 



-1 



e A 

-A* -e 



(17) 



where S is an infinitesimally small positive number. Using 
this representation for g^ it is easy to derive boundary 
conditions for Eq. (|l5|). In a normal reservoir A = 
and g^ — — o z . In a superconducting terminal having 
phase (j), A = 

•A 



i<t> 



and the Green function satisfies 



|A S |e- 
g™ = <r x sin (p + dy cos <p. 

It is also possible to simplify Eq. ( fl3"l) for the distri- 
bution functions considerably: In the case of a negligible 
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supercurrent I s = J j e de the equations for the two distri- 
bution functions decouple, reducing /l to its equilibrium 
value /l = tanh(^) and leaving us with a single equa- 
tion for /», which can be cast into the form of a diffusion 



equation:li3 

V (D(e, r)V/ T (e, r)) - 7 (e, r)/ T (e, r) = 0, (18) 

where the first term describes diffusion of quasiparticles 
with an effective diffusion coefficient 



V 

1 

V 



Tr{i - a z g?a z g*}, 



= ^-Tr{& A + <^) 2 }, 



(19) 



that is modified by the penetrating superconductivity. 
Here we have used the identity gf- = —a z g^ £ a z , which 
relates the advanced and retarded Green functions. The 
second term describes absorption of quasiparticles in 
the superconducting condensate with a coefficient 7 = 
iTr{(gf + <7 A )A}. In the absence of external fields 7 is 
proportional to the local value of the pair potential, since 
in that case we can always choose A to be real, and we 
obtain 



; +9?)}- (20) 
Jl8| ) follow from ex- 



7 (r, £ ) = --A(r)Tr{^ 

The boundary conditions for Eq. 
panding tanh( ) to first order in V. This deter- 
mines the boundary condition for /t . In a normal reser- 
voir that is biased at a (small) voltage V with respect 
to a superconducting lead, the distribution function is 

./T = (#)cOsh- 2 (^). 

In most theoretical approaches, see e.g. Refs. [ 
||,|],[lo]|26) , electron-electron interactions in the normal 
metal are disregarded, leading to A, 7 = 0. However, 
as shown in Ref. [l6|, including the effect of these interac- 
tions produces a change in the resistance. We thus need 
Eq. ( |l4| ) for A, which after a straightforward calculation 
can be rewritten as: 



A = 



A 



de tanh 1 



— 1 Tr {ia y (gi 



277 



gf)}- (21) 



We now have all the necessary ingredients to calculate 
the various non-equilibrium transport properties of the 
system. The next section will be devoted to two of these 
properties, namely the electrostatic and chemical poten- 
tial distributions. 



IV. 



ELECTROSTATIC VERSUS CHEMICAL 
POTENTIAL 



As a first example of the theory of section III we cal- 
culate the electrostatic and the chemical potential in the 
ID wire of Fig. ||(a). As was shown above, the former is 
determined by the distribution of electric charge in the 
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c) 

FIG. 2. (a) The simplest possible system, consisting of a 
diffusive normal metal wire of length L, that is connected 
on the left to a normal reservoir and to the right to a su- 
perconducting one. (b) An example of a more complicated 
structure. The branches I, II and IV have length L', branch 
III has length L. For further details see text, (c) The coor- 
dinates on the branches. The origins are indicated for each 
branch. 



wire. The latter determines the magnitude of the cur- 
rent that flows through the sample. The system consists 
of a diffusive normal metal wire of length L attached on 
the left to a normal metal reservoir and on the right to 
a superconducting terminal. The normal lead is biased 
at a small voltage V with respect to the superconductor. 
The normal reservoir is situated at x = 0. 

Using Eq. ( |l2| ) for the Kcldysh component and Eq. ( |T6| ) 
for the Green function we obtain 



1 



de /t(x,£) cos(0(x, e)). 



(22) 



In Fig. y we have calculated the potential distribution 
in the normal metal for different values of L/£, i.e. for 
different temperatures. 

In the limiting cases of low and high temperatures, the 
potential can be calculated analytically: 



<p{x) 
ip(x) - 



.V[l 

v(i- 



X\ [1TX 

— cos — 
LJ \2L 



V 



< k B T < A s 



(23a) 
(23b) 



Fig. H shows that the potential distribution changes with 
temperature from a non trivial one which is influenced 
strongly by the penetrating superconductivity to the ex- 
pected linear dependence for high temperatures (but still 
k^T -C As). This behavior is caused by the fact that the 
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FIG. 3. Electrical potential distribution in the wire of 
Fig. |^(a) as a function of temperature. Going from bottom 
to top, the curves correspond to values of L/£ = 0, 1.0, 1.5, 
2.0, 4.0, and oo. The temperature is proportional to (L/^) 2 . 

density of states vanishes in the vicinity of the supercon- 
ductor. Hence, the charge distribution which causes the 
electrostatic potential also vanishes in this region. The 
most important consequence of the nonlinear voltage dis- 
tribution across the sample is the fact that the resistance 
at a certain point is no longer local, but depends on the 
distribution of resistivity in the entire structure. This is 
a direct consequence of the coherent nature of Andreev 
reflection as was discussed above. 

The chemical potential, which is simply proportional to 
the energy integrated distribution function, is much less 
sensitive to changes in temperature. The zero and high 
temperature distributions are the same and are given by 

l*o(x) = v(l-j). (24) 

In Fig. U we have plotted the deviation of the chemical 
potential from this zero temperature solution. First of 
all we note that the change is very small. The maximum 
change at first increases rapidly with increasing temper- 
ature. However, beyond a certain temperature, L/£ » 4, 
the maximum starts decreasing slowly. For high tem- 
peratures, the solution returns to the linear distribution. 
Thus we see that the electrostatic and chemical poten- 
tial are only equal in the case of high temperatures, when 
both reduce to the trivial linear dependence that is also 
exhibited in the case of two normal terminals. 

To indicate more clearly the difference between elec- 
trostatic and chemical potential, we have plotted the ra- 
tio of the two for different temperatures in Fig. §. The 
electrostatic and chemical potential differ most at low 
temperatures and near the superconductor. However, in 
the vicinity of the normal reservoir they are almost equal 
and the extent into the wire at which they are equal in- 
creases at higher temperatures. In the high temperature 
limit they both reduce to the same linear distribution, as 
was shown above. 



0.00 



-0.03 



1 1 1 1 1 1 1 1 1 












\\\\ 1.5 / 








V\\\\ 2.0/ 










4.0 / 










/ 6.0 




i,i, 


, 50 , 


i 



0.0 0.2 0.4 0.6 0.8 1.0 

x/L 



FIG. 4. Change in the chemical potential for the wire of 
Fig. ^(a) relative to the (linear) zero temperature distribu- 
tion as a function of temperature. The labels near the curves 
indicate the value of L/£. 




FIG. 5. The ratio of electrostatic and chemical potential 
as a function of temperature. Going from bottom to top, the 
curves correspond to values of L/£ = 0, 1.0, 1.5, 2.0, and 4.0. 

V. THE INTERACTION EFFECT 
A. The resistance of a ID wire 

In this section we calculate the resistance change at 
zero temperature due to the penetration of the pair po- 
tential into the normal metal region. We first consider 
the wire of Fig. ^(a) and then address more general ge- 
ometries. In order to simplify the equations later on, the 
superconducting end of the wire is now located at x = 
and the normal end at x — L. 

We first solve Eq. ( |l5| ) for g A numerically to zeroth 
order in A and ip : 

PV (g A V.g A ) - ie[a z , ,g A ] = 0, (25) 

where ip = (zero voltage limit) and A has been disre- 
garded since it is small in the normal metal region be- 
cause of its proportionality to the interaction parameter 



G 



A. We then solve Eq. ([l8]) to first order in 7. At zero tem- 
perature, D(s, x) = T>, and we write /t = /o + fi with 
/1 < f . Then f = c(e)§ and /{' = ^/ , with c(e) = 

|^cosh _2 (^). Using the fact that J Q L f[(x')dx' = 
(A(0)=/i(£) = 0) we find 



/o(i) 



£(£) 
L ' 

£(£)_ 

VL 2 



x j(x)dx. 



(26a) 
(26b) 



The current in a normal piece of metal is proportional 
to the local gradient of /. Since everywhere in the dif- 
fusive wire there is an induced pair potential it is not 
possible to calculate the current in this way somewhere 
in the middle of the wire. We must calculate it at the 
normal reservoir, where the pair potential is forced to 
vanish. The current flowing out of the normal reservoir 
(and hence the conductance G of the system) is therefore 
proportional to f'(L). Noting that = — ^ we obtain 
for the relative resistance change: 



-R- = -7m = -vLj dxx 1{xl 

where 7 is given by 

7(2;) = 2A(x) sin(0(e = 0,xj). 



(27) 



(28) 



Thus the relative resistance change is proportional to the 
interaction parameter A and its sign depends on the sign 
of A. Here we use the convention that A is positive for at- 
tractive effective interactions in the metal. Furthermore, 
the resistance change depends sensitively on the precise 
geometry of the structure, as will be shown in more de- 
tail in the next section. A measurement of this resistance 
change would allow one to directly measure A. Calculat- 
ing the relative resistance change for the wire of Fig. ||(a) 
gives 



(29) 



independent of T> and L. For silver, the estimated value 
of the interaction parameter is A = +0.04, (Ref. [2^ ), and 
therefore the resistance of a silver wire in contact with 
a superconductor is reduced by 5.5% with respect to its 
normal state value. 



to the one we used recently0 to model the experimental 
setup of Ref. [l4|. In the experiment a current was applied 
along the path indicated by the arrow going from normal 
lead N\ to normal lead N 2 - The voltage was measured 
between the opposite normal reservoirs A3 and A4. The 
superconductors have a phase difference 4>. In principle 
Eqs. (JT^) and ([l8]) should be solved in each branch of the 
structure and the solutions matched in every nodal point. 
However, it is possible to reduce this geometry to a one 
dimensional one due to the symmetry of the geometry: 
The voltage and /t distributions are antisymmetric with 
respect to the line Si — S2 in Fig. ||(b) whereas A and 
the Green function g e are symmetric. This allows us to 
consider only the three elementary branches I, III, and 
IV in the calculation of g e . 

For the calculation of the resistance change we need 
only consider the branches I and II, each having a length 
L' , and III which has a length L. As depicted in Fig. 0(c) 
we use coordinates x, y and z to denote the position on 
branch I, II and III, respectively, with the origins as 
indicated. The zeroth order solution of Eq. (jig) is 




fa 



Integrating /{' = ^/q in each branch gives 



(30) 



ci + /„ dx' (L + L' 
fi = { c 2 + j y dy'L2m 
c 3 + f*dz> z'^f 



x') 



7(3=') 
v 



(I) 

(II) 

(HI) 



(31) 



Here the integration constant ci = because the cur- 
rent in point x = is fixed and c 2 — because there is 
no current in branch II. Note th at th e situation here is 
different from the one in section VA. There we calcu- 



lated the change in current at fixed voltage whereas in 
this case we calculate the voltage change at fixed cur- 
rent. Consequently the resistance change is now given 
by ^ = ^jr — jH^ . In the nodal point A we have 
f[(x = L 1 ) +f[{y = L') +f[( z = L)=0 and this gives 

1 f L , *7(«) , f L \ Lj(y) 

C3+ L dz — + L dy — 



V 



(32) 



B. Generalization to arbitrary geometries 

The results obtained above are readily generalized for 
arbitrary structures containing a number of normal reser- 
voirs, two superconducting ones and a number of diffusive 
wires connecting them. To illustrate how the resistance 
change for such systems is calculated we consider the ge- 
ometry shown in Fig. ||(b). This is a structure similar 



The distribution function /1 in point A is now given by 
h(A)=c 3 [ L dz+ [ L dz f dz' Z ' 1{Z '^ 



c 3 L 



dz 



Jo 
Lz - z< 



V 



V 



7 (z). 



Analogously we find for fi(B): 



(33) 



7 



h{B)=h{A)- [ dy I ,lu 
Jo Jo 



v 



L y 2 - LL' 

dy — ^ — i{v)- 



V 



(34) 



Usin g E q. ( [32] ) to eliminate the first two terms in expres- 
sion (|3^) for fi(A) we calculate the resistance change: 



6R 
~R 



V 



L-U 



v 



dz 



z 2 l{z) 



VL 



dx 



(L + L' - x)j(x) 
V 



(35) 



The relative resistance change depends on the phase dif- 
ference <p between the superconductors through 7, which 
in its turn depends on <fi via the boundary conditions 
imposed on g E . Note that although branch IV does not 
enter the calculation of the resistance change, its length 
does influence the Green function, hence 7, and thus also 
the resistance of the structure. 

Before discussing the results of a calculation for an ex- 
perimentally relevant geometry, it is instructive to look 
at the (slightly unrealistic) case of constant 7 to investi- 
gate the qualitative behavior of the system. In this case 
the resistance change reduces to 



SR 



7^ 
V 



1 3 
-a J 

3 



3 1 \ L' 

—a 2 — 2a I with a = — (36) 

2 3 / L 



which shows that, since 7 ~ A-p-, the relative resistance 
change is proportional to the interaction parameter and 
independent of the absolute values of L, L 1 and T>. Only 
the ratio of L and L' is important. In this simple case 
the largest effect would be obtained for a — 3.56. We can 
get an idea how sensitive the resistance change is to the 
details of the geometry by comparing the cross structure 
{L 1 = 0) with the structure in which V = L. It is easily 
calculated that the resistance change is higher in the lat- 
ter case by a factor of 10.5, showing that a relatively small 
change in geometry causes a major modification of the 
resistance change. In a realistic computation the modifi- 
cation due to these side branches is not so dramatic, but 
is still about a factor 3. 

The result of such a realistic calculation is depicted 
in Fig. ^ where the scaled resistance change has been 
plotted for the Andreev interferometer of Fig. ^(b) with 
Li = Li 1 = Liv = L and Lm = 2L. This is the layout 
we used in Ref. 16 to model the experimental setup of 



Ref. _ 
Petrashov et al. 

There are two main differences with the previous re- 
sults for the one dimensional wire. The first is the depen- 
dence of the resistance change on the phase difference <j> 
between the superconductors with oscillation period 2n. 
The effect has its largest (negative) value for zero phase 
difference and vanishes for a phase difference of 7r, when 
superconductivity in the current branch (branches I, II, 
and their mirror images) is completely suppressed. The 




FIG. 6. The scaled resistance change due to the in- 
teraction effect of the structure shown in Fig. ^(b) with 
Li = La — Liv = L and Lm — 2L. The magnitude of 
the effect is 0.6%. 



second difference is the magnitude of the effect. In the 
one dimensional case this was 5.5%. Here it is 0.6% (for 
silver). Although this value is proportional to the inter- 
action parameter and strongly depends on the length of 
the branches I and II, as was established above, a more 
important reason in this case is that the superconductors 
are not in the current path. As a consequence, the pair 
potential in branches I and II is smaller than in the one 
dimensional case, leading to a weaker effect. In general, 
the magnitude of the interaction effect is smaller than the 
influence of a finite temperature on the resistance. This 
latter phenomenon is the subject of the next section. 



VI. THE THERMAL EFFECT 

A. The temperature-dependent resistance of a ID 
wire 

As was shown recently by the authors,0 the exper- 
imental results of Ref. [l4| could be explained by the 
temperature-dependent proximity effect. In this case no 
pair potential is induced in the normal metal region by 
the electron-electron interaction but coher ence occurs as 
a result of the finite extent £ = ^JV/T at which the 
superconductivity penetrates into the normal metal. Al- 
though the implications of this effect for complex struc- 
tures are not always immediately apparent, the mech- 
anism itself is net new. It was already studied in the 
middle seventiesE-3 and is in fact the phenomenon of An- 
dreev reflection in diffusive metals. In spite of the fact 
that the effect has been known for a long time, a clear 
physical picture is still lacking. 

To describe this effect we disregard the second term 
in Eq. (pit) but we now take into account the fact that 
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1.00 



at finite temperatures the effective diffusion coefficient 
is no longer constant but depends on energy and posi- 
tion: D{e,x) = fTr{({£(x)+gtt(x)) 2 }. The diffusion 
coefficient reduces to T> for low and high energies. For 
energies e f=a D/L 2 it exhibits a maximum which is about 
twice as large as the zero energy value. The temperature 
enters the boundary conditions for Eq. fll8|), but more 
importantly it determines the energy window in which 
the quasiparticles experience the energy dependence of 
the diffusion coefficient. Therefore the resistance change 
should vanish at both low and high temperatures and to 
have a minimum for quasiparticle energies e ss D/L 2 . 

As in the previous section, we will first calculate this 
effect for the one dimensional structure of Fig. ^(a) and 
then extend the treatment to more general structures, us- 
ing the geometry of Fig. ||(b) as an example. The normal 
lead is now situated at x = and the superconducting 
one at x = L. The procedure is as follows: We first in- 
tegrate the diffusion equation for /t twice, which gives 
/t(x, £•) = a(e) f Q D~ 1 (x',e)dx' + 6(e). The integration 
constants may still depend on e and are determined by 
the boundary conditions /t(0, e) = |yCOsh _2 (^) and 
h(L,e) = 0. This gives 

h(X) £) = z_V 2 / e \ / _ rn^e) \ 

3 v ' ' 2T \2TJ V m(L,e) J v ; 



with 



m(x,e) = — / D 1 (x',e)dx'. 



(38) 



As before, the current flowing out of the normal contact 
is proportional to /^,(0,e) and hence we obtain for the 
normalized conductance 



G 

G N ' 2VT 

e 



'If 00 / E \ 

—— I de m~ 1 (L, e) cosh -2 ( — - ) 

de m -1 (L,e)cosh -2 (jfif) , ( 39 ) 



2L 2 V 



where we have used T = T>/^ 2 and scaled the energy 
with the Thouless energy e = eL 2 /T>. Note that the 
cosh~ 2 (^y) part of the integrand defines the energy in- 
terval in which the quasiparticles feel the energy depen- 
dence of m -1 (L,e), as mentioned above. Eq. ( p9[ ) shows 
that the normalized conductance depends on the ratio 
L/£ ~ \/T only (the factor [Dm(L, e)] _1 is dimension- 
less and does not depend on temperature). In Fig. |?] the 
temperature dependence of the normalized resistance, 
= ^j-, is plotted. As expected, the effect vanishes 
at T — and for T — > oo. The maximum magnitude 
of the temperature effect is material independent and in 
general larger than that of the interaction effect; 10% vs. 
5.5% for the particular case of a silver wire. Although 
the high temperature tail like the one of Fig. |t] was al- 
ready observed in Ref. the low temperature reentrant 
behavior of the resistance has not yet been measured. 



0.98 



0.94 



0.92 



0.90 



10 
L/c 



FIG. 7. Temperature dependence of the normalized resis- 
tance of the structure shown in Fig. ^(a). The temperature 
is proportional to (L/£) 2 . 



B. Temperature effect in Andreev interferometers 

A generalization of this theory to more complicated 
structures should take into account how a resistance mea- 
surement is actually done. In the case of the experiment 
of Ref. [l4| a current was applied between N\ and N2 
whereupon the voltage between A^ 3 and was measured 
(see Fig. As a consequence, there is no current in 
branch II nor in its mirror image attached to reservoir JV4. 
To calculate the resistance of this complex structure we 
proceed as follows: We start by noting that the current 
is conserved for each quasiparticle energy e. This implies 
that Ti{e) +Iu{e) — liu(e). We now apply standard 
circuit theory for every quasiparticle energy (we use the 
fact that due to the antisymmetry of the voltage distri- 
bution the voltage is zero in the middle of the structure): 



Vi(e) =li(e)Ki(s) +l IU {e)n ni {e), 
V 3 (e) =ln(£)1l u (e) +l IU (e)Ki U (e), 



(40a) 
(40b) 



where V,(e) = ^ cosh - (^) and the indices on V refer 
to the respective reservoirs and the indices on X and 1Z 
refer to the appropriate branches. The proportionality 
constants TZi.iuiiis) in the different branches can be 
deduced from the 'resistance' of the branches connecting 
two reservoirs: 



Hi4 (e) 



f'M-f'M 



(41) 



where i and j label the reservoirs at which the distribu- 
tion functions are evaluated. 

Using the current conservation condition, Eq. ( |40| ) 
can be cast into the form of a matrix equation V(e) = 
lZ(e)l(e), where TZ(e) is a (2 x 2) matrix relating the 
voltage and current vectors V(e) = (Vi(e), V3(e)) T and 
1(e) — (li(e),Iu(e)) T . Inverting the matrix equation 
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FIG. 8. Scaled resistance change due to the thermal effect 
of the structure shown in Fig. ^(b) with Li = Ln = Liv = L 
and Liu = 2L for L/£ = 3. The magnitude of the effect is 
9.7%. The squares are the experimental data of Ref. [l4|. 



and integrating over all energies gives / = GV, where 
V = (Vi,V3) T , / = (Ii,Iu) T , and the conductance ma- 
trix is given by G — ^ deTZ~ 1 (e) cosh~ 2 (^). We 
now impose the boundary condition mentioned above on 
the currents in branches I and II: I = (I, 0) T . This gives 
a relation between V\ and V3 in terms of the matrix ele- 
ments of G: V3 = —tt^Vi, which in turn must be used to 
calculate the experimentally measured resistance of the 
structure: R = 2ji = 2(G U - G g^ )~ 1 . The resis- 
tance depends again on the phase difference between the 
two superconducting reservoirs because of the boundary 
conditions imposed on the Green function. 

In Fig. H we have plotted the results of this calcula- 
tion for the previously used geometry of Fig. ^b) with 
Li = Li 1 = Liv = L and Lm — 2L. This particu- 
lar calculation was done for the value of L/£ = 3 which 
gave the biggest effect. Also shown are the experimental 
results of Ref. |l4|. 

Again we find 2n periodic oscillations like in Fig. [fj, 
but although the phase dependence looks similar there 
are some distinct differences. First of all the shape of 
the oscillations is different, especially near the minima, 
where the thermal effect produces less narrow peaks than 
the interaction effect. The amplitude of the thermal os- 
cillations is material independent whereas that of Fig. ^ 
is proportional to the interaction parameter of the metal. 
Moreover, the maximal magnitude of the oscillations, 
9.7% in this particular case, are in general much larger in 
the thermal effect, which makes observation of the pre- 
viously discussed oscillations more difficult. 

We obtain an excellent fit with the experimental data 
in Fig. ||, where the magnitude of the oscillations is about 
11%. However, a decisive check on our theory would be 
provided by the observation of the temperature depen- 
dence of the resistance like the one plotted Fig. fj] for the 



system considered here. Indeed, unpublished data by the. 
authors of Ref. [l4| also, show a remarkable agreements 
with our calculations .ta However, the reentrant behavior 
has not yet been observed unambiguously. 



VII. SUMMARY 

We have performed calculations of the distribution of 
the electrostatic and nonequilibrium chemical potential 
in a one dimensional disordered superconducting hybrid 
wire. We showed that the two behave differently as a 
function of temperature and that they are in fact only 
equal in the high temperature limit. We have proposed 
an experimental setup to measure these different poten- 
tial distributions. We have also presented-a-computation 
for this wire using a recently developedliS mechanism, 
that causes the resistance of a diffusive superconducting 
hybrid structure to change at zero temperature. The lat- 
ter is in contrast with the well known thermal mechanism 
of Andreev reflection in diffusive metals, where the resis- 
tance change vanishes for low temperatures. In addition 
to this new result we have given a detailed account of 
the calculation performed in Ref. |l6| for the experimental 
setup of Ref. [l4 . Since the relative resistance change due 
to the novel mechanism is proportional to the interaction 
parameter in the normal metal, observation of this effect 
would allow a direct measurement of this physical quan- 
tity. Furthermore, we have shown how to calculate the 
resistance of an arbitrary Andreev interferometer using 
finite-temperature proximity effect theory. Because the 
finite temperature effect generally causes a much larger 
resistance change, this is the correct theory to describe 
an experiment like the one performed in Ref. 14. 
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